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The self gravity of an infalling gas can alter significantly the accretion of gases. In the case of spherically 
symmetric steady flows of polytropic perfect fluids the mass accretion rate achieves maximal value when the 
mass of the fluid is 1/3 of the total mass. There are two weakly accreting regimes, one over-abundant and the 
other poor in fluid content. The analysis within the newtonian gravity suggests that selfgravitating fluids can be 
unstable, in contrast to the accretion of test fluids. 



INTRODUCTION. 

The analysis of the spherical steady accretion of gas onto 
a gravitational center starts with the Newtonian accretion dis- 
cussed by Bondi Michel |2], Shapiro and Teukolsky 
and others jGI] investigated the accretion of test fluids in the 
Schwarzschild space-time. Recent analysis showing the im- 
portance of backreaction is that of Karkowski et al. 1 5 ] . In this 
lecture we follow the general-relativistic description of self- 
gravitating fluids that has been formulated in |6] and present 
main results of i5|j] . 

One of the reasons for inspecting steady flows in the full 
generality, when infalling gas modifies the space-time geom- 
etry, is the wish to see the effects of backreaction in a simple 
but nontrivial accretion model. We show later, following 
that the mass accretion rate m behaves like (1 - xf-niB, where 
X - nif/m is the ratio of the mass of the fluid to the total 
mass and riiB is the mass accretion rate for the accretion of 
test fluids. Since ihs is proportional to the fluid abundance, 
hence to x, one sees that there exists a maximum of ih when 
X - 1/3. The mass accretion rate m sets an upper limit for the 
luminosity of the system. This simple analysis of accretion 
suggests the existence of two weakly luminous regimes: one 
rich in fluid with nij ^ m and the other with a small amount 
of fluid, nif/m <K 1. It is clear from this reasoning that it 
is reasonable to take into account selfgravitation of accreting 
gases whenever their mass is comparable to the central mass. 
On the other hand, the stability properties of accreting fluids 
can depend on backreaction, however small. The stability of 
steady flows in a fixed relativistic background as well as in 
the newtonian context has been shown by Moncrief |7]. We 
show later that the steady newtonian acretion of selfgravitat- 
ing gases does not exclude instability. 

Trying to put matter in the astrophysical context, notice that 
one can expect the mass of accreting fluids to be comparable 
to the central mass only in the case of heavy and luminous 
compact objects; compact galactic centers probably constitute 
fair examples. This suggests that the accretion physics of such 
objects should include the selfgravity of accreting matter Let 
us mention here that there exists a fundamental physical as- 
pect that the astrophysics of accreting matter might be able 
to resolve. Namely, the standard point of view is that com- 
pact objects of a mass exceeding several solar masses consti- 
tute black holes. That interpretation has been challenged by 
Mazur and Mottola f3] who developed the notion of gravatars 



— stars almost as compact as black holes but dispossessed 
of event horizons. Gravastars are supposed to be built from 
matter violating energy conditions |9] and in this way they 
would avoid the standard compactness and mass limitations, 
like the Buchdahl theorem or Chandrasekhar type bounds 
on mass. There is an interesting motivation coming from the 
quantum-gravity type arguments, that inspired the concept of 
gravastars, but we leave it aside. It is disputable whether ob- 
servations can distinguish compact bodies endowed with a 
hard surface from black holes within the known nowadays ac- 
cretion models. For a discussion of that issue see |11] and 
We believe that the resolution of that problem requires 
the inclusion of general-relativistic aspects such as backreac- 
tion. In this talk we shall report the recent progress in this 
direction. 



MAIN EQUATIONS. 

We will consider a spherically symmetric compact ball of a 
fluid falling onto a non-rotating black hole. The equations are 
taken from We will use comoving coordinates 



ds^ = -N^dt^ + adr^ + [dO'- + sin^ Od^'^) , 



(1) 



where the lapse A^, a and the areal radius R are functions 
of a coordinate radius r and an asymptotic time variable t. 
The nonzero components of the extrinsic curvature Kij of the 
t = const slices read K' = ^^5,0-, = Kf, = l§ = 

jitiK - /T^). Here tiK - jjdt In ^ -s/aR^y The mean curvature 
of two-spheres of constant radius r, embedded in a Cauchy 

2(5 R 

hypersurface is A: = The energy-momentum tensor of 

the perfect fluid reads T'"' = {p + Q)d^u' + pg^^'^ (where we 
assume the barotropic equation of state), m^m'^ = -1 with m'' 
denoting the four-velocity of the fluid, p is the pressure and g 
the energy density in the comoving frame. 

The areal velocity f/ of a comoving particle labelled by co- 
ordinates (r, f) is given by U{r,t) - - j(\xK - K'^). From 
the Einstein constraint equations one has 



kR^2 



R 



(2) 



Here and below we use dr - ■\/a(kR/2)dR in order to elim- 
inate the differentiation with respect to the comoving ra- 
dius r. The quasilocal mass m(R) is defined by dRm(R) = 



2 



AnR^g. The mass accretion rate is m - (5, - (d,R)dR)m(R) - 
-4-nNR^U(g + p) or (equivalently, but perhaps more famil- 
iarly) m - -AnR^nU, where n is the baryonic density. We 
assume a steady collapse of the fluid, which means that 
all its characteristics are constant at a fixed areal radius R: 
5,X|R=co„st = {d, - (d,R)dR)X = 0, where X = g,U,a^. 
a = ^J^pp is the speed of sound. Strictly saying, a station- 
ary accretion must lead to to the increase of the central mass 
and of some geometric quantities. This in turn means that the 
notion of the "steady accretion" is approximate — it demands 
that the mass accretion rate is small and the time scale is short, 
so that the quasilocal mass m{R) does not change significantly. 
Thence also geometric quantities, like the mean curvature k or 
the area of the black hole, remain practically constant. 

The Einstein evolution equation drU = k^R^duN/A- - 
m(R)N/R^ - AnNRp and the energy conservation equation 
dtQ - -NtiK{g + p) constitute (assuming steady flow) ordi- 
nary differential equations with respect to R. 

It is easy to show that now the mass accretion rate is con- 
stant, dRfh -Q^. 

By a black hole is understood a region within an apparent 
horizon to the future, i.e., a region enclosed by an outermost 
sphere 5 a on which the optical scalar 6+ = ^ + U vanishes 



c» = InRlp^. The infall velocity U reads 



(The other condition, that 0_ = ^-U > for all spheres 
outside Sa, is satisfied trivially for a steady accretion.) That 
means that on Sa the ratio ImbhIRah becomes 1, where Rah 
is the areal radius of the apparent horizon and iiibh = '«(^ah) 
is the mass of the black hole. We assume the polytropic equa- 
tion of state p - Kg^, with F being a constant, (1 < F < 5/3). 

Furthermore, we assume that the radius Rcc of the ball of 
fluid and other boundary data are such |f/co| <K 'Jl^^ <g 
Qoo. These boundary conditions allow one to join smoothly 
the internal solution (the steady fluid) with the external 
Schwarzschild geometry. The momentum conservation equa- 
= — in comoving coordinates 



tion V^r^^ 



NSrp + ip +g) OrN = 
— can be integrated, yielding (with N{Roo) =1) 



-F + 



r + al 



N" 



(3) 
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Here k - Equation (0} is the general -relativistic version 
of the Bernoulli equation. Equation Q implies also that the 
lapse expresses as follows - with the baryonic num- 
ber density n = Cexp^'^ ds-^-^^^. The baryonic current is 
conserved, (nu'') = 0. 

The sonic point is as a location where the length of the 
spatial velocity vector \U\ = \U\/{kR/2) equals a. Therefore 
\U\ - jkRa. In the Newtonian limit this coincides with the 
standard requirement \U\ = a. In the following we will de- 
note by the asterisk all values referring to the sonic points, 
e.g. a,, Ut, etc. The four characteristics, a,, t/*, ot,//?, and 

are related H, flj (l - ^ + c.) = f/2 = ^ -h c., where 



U 



1 -1- 



(5) 



Here f/» is the negative square root. From the relation be- 
tween the pressure and the energy density, one obtains, using 
equation (0} 



g - gca 



-I) 



F 



^ + 1 



(6) 



where the constant g^o is equal to the mass density of a col- 
lapsing fluid at the boundary Roo- The steady fluid is described 
by equations (|4]-E}. They constitute an integro-algebraic sys- 
tem of equations, with a bifurcation point at the sonic point, 
where two branches (identified as an accretion or a wind) do 
cross; that is a well known feature of that problem, present 
also in models with a test fluid CI - 01 > El- That requires 
some caution in doing numerics and a careful selection of the 
solution branch. Notice that these equations are expressed ex- 
clusively in terms of quantities that are steady in the sense of 
the former definition. 



CALCULATION OF THE MASS ACCRETION RATE. 

In numerical calculations it will be convenient to represent 
the mass m{R) in the form 



m{R) 



— m - An \ I 

JR 



drP'g. 



(7) 



Notice that ;«(oo) - nice ~ m. For a steady flow one obtains 
L6J 



kR 

N = .„ /3iR), 



kiR^)R. 
/3(R) = exp(-16;r 



Jr ^'"-'^¥1 



(8) 



The full description of the steady accretion is given by equa- 
tions 0, (0}, and (jSJl, which are equivalent to the for- 
mer set of equations. 

In the forthcoming part of the talk we shall describe main 
points of the derivation of the accretion mass rate, following 
the presentation of Karkowski et aljpl]. 

In the first step, one shows in pi], that outside the sonic 
horizon solutions of the equations with backreaction can be 
estimated by relevant solutions describing the stationary ac- 
cretion of test fluids. We supplement the foregoing list of con- 
ditions by another, that the mass ntf of the fluid outside the 
sonic horizon is of the order of the mass m, within the sonic 
horizon; That allows one to show that the function yS appearing 
in equation (|8} is essentially 1 and c, = InR^p^ <Si al can be 
ignored. As a consequence one infers that in the case of sta- 
tionary accretion of selfgravitating fluids some characteristics 
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of the sonic point — a J, Ul and m^/R, — are the same as for 
the test fluid accretion, provided that asymptotic data p^o and 
are identical |5]. In particular, in both cases characteristics 
of the sonic point do not depend on the mass — either the total 
mass m or the mass m» of the accretor 

In the second step one shows that it is possible to estimate 
the energy density p of a system with backreaction using in- 
formation about solutions in the test fluid approximation. It 
is done in 01 . We shall describe in a more detail how one 
obtains the desired form of estimates onto the energy den- 
sity of test fluids. One starts with the estimate on the inverse 
of the lapse function (6) - VT^, where 5 s ^ - [/^ 
is less than one (and in fact strictly smaller than 3/4). It is 
easy to show, applying the method of differential inequalities, 
that N-^ (6) < I + SI {iN^ {6)). Since 5 < 2m/R, one has 

N-^ (6) < I + ml [rN^ (5)). Notice now that the lapse is an 
increasing function and N(R) > Nt = N(R») for R > R^. 
Therefore one arrives at the inequality 



1 m 
< 1 -H — X -. 

N^. R 



(9) 



Having that, one can derive a bound onto g. One obtains, 
applying the inequality (|9j. 



eiR) < 



1 + ^ 



i-f(r-i)- 



R 



i/(r-i) 



(10) 



The lapse at the sonic point A^* = 1/ -\/rT"3flf and m/7?» - 
oiNl; thus a constant a = 



fn: r. 



can be written as 



iJl+3alial + ^ 
1 \ aZo 



One can prove following bounds (see equation (5.2) in 

^2 2 



ai, 5 . 



•3rH- 



3fl;(r-iy(9r-7) 

4r(T+3fl!) 



(11) 
1) 

(12) 



This implies a series of estimates: i) for 1 < F < 4/3 one 
has < 2; ii) in the range 4/3 < T < 3/2 holds ^ < 4; 

iii) when F > 3/2 and < 0.2 one obtains ^ < 2.3/aJ and 
flj < 1.6flco, iv) while for F > 3/2 and > 0.2 it follows 
^ < 4. Now, if we demand that the sonic point is outside of 

an apparent horizon then m/ (2/?») < 1 /4 and that implies a J < 
1 . Using this fact, one can derive bounds onto the constant a. 
Namely, in each of the above cases we have i) 1 < F < 4/3 — » 
a < 6; ii) 4/3 < F < 3/2 ^ a < 9.5; iii) F > 3/2 and 
al < 0.2: a<\+ 2/a<» iv) F > 3/2 and a] > 0.2: a < 8.3. 

With these estimates one obtains simple bounds onto p. We 
shall consider in more detail two of them corresponding to the 
cases i) and iii). For F < 4/3 equation J10> can be written, 
using the inequality {\ + AIxY < and a < 6, as 



q{R) < Qoc exp 



6R, 



(13) 



1 

In the case iii) one notices that (1 -I- A) r-i < 1+2A+A . That 
leads to the estimate 

P R* 
— < 1 -h2(F- l)^— -1- 

Poo R 



(r-.,'»'|. 



(14) 



one can check the selfconsistency of the calculation — for 
instance that the function /3(R) * 1 for R > R,. Indeed, the 
exponent in the exponential funtion yS (Rt) is of the order of 
Mf = An drrp. This latter expressions reads in the case i) 



Mf 



2np^R 
6R, 

Rex, 

36Rl 

rL 



Ri 



Rm 
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(15) 



Since R^a » /?» one concludes that Mf » 2npcoR\a which 
in turn is clearly much smaller than mIR,. Thus the expo- 
nent of f3(R) is negligibly small for R > R^ and li(R) ^ 1 
as desired. A similar reasoning in the case iii) leads to 
Mf ^ 27Tp^Ri(l+0(aRJRo,)). Now aRJR^ < RJR^ + 
2Rtl (Rooaca). It is easy to see, taking into account the bound- 
ary condition » m/7?oo and a, ^ Ak,, that the second term 
in Mf is much smaller than 1 and therefore the expression 
O (aR^IRao) is negligible small. Thus again Mf a; niflRoo <^ 
mIRca <K 1, where irif is the mass of fluid outside the sonic 
horizon; therefore one again arrives at /3 (/?,) ~ 1. 

In the third step one shows — using the same reasoning as 
in the preceding paragraghs — that the mass nif - -/Qoo, where 
y is a constant. This is because the mass density changes mod- 
erately outside the sonic point, as shown by estimates ( I13ll4t . 
The calculation is lengthy but simple and we omit it here. Tak- 
ing into account m - m^+nif, one infers m,/m - l-Qoaiylm); 
the ratio x = m, / m linearly decreases with the growth of ■ 

The rate of the mass accretion rii can be expressed as below 
(see equation (6.1) in |6]) 



2 / « 
tn - mn^Qm — — 

mi \ai 



2 \ 2(r-i) 



1 -I- 



1 +3ai 



(16) 



The whole dependence on g^o is contained in the factor 
mlgoo - m^xi l - x)^ ly. Thus one can write 

(5-3r) 

,n^mx——n-J-^ 1+ (17) 

This expression clearly demonstrates that the mass accretion 
rate achieves a maximum at m* = 2m/ 3 and tends to zero 
when i) ;«» — > m (when the density g tends to zero) and ii) 
m^tlm — > 0. 
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Thus, one concludes that Amongst steadOy accreting sys- 
tems of the same Rc^, Ooo and m those will be most efficient 
for which to, = 2to/3. The factor 2/3 is universal — indepen- 
dent of the parameters 7?oo, T and aoo- 

Since the mass of the fluid to / is close to to - to, and to - 
nif + niBH, the above means that the maximum of the mass 
accretion takes place when the mass of the fluid is a half of 
the mass of the black hole. 

Finally let us mention that in the case of polytropes with - 
Krf the Bernoulli equation reads (t - 1 - a^) = F - 1 - a^. 
The analytic arguments of the type presented above can be 
applied and conclusions are the same. Similar results should 
be true for massive fluids in the newtonian limit provided that 
the poly tropic index F is clearly smaller than 5/3. 

NUMERICAL RESULTS. 

In this Section we present a small sample of the numerical 
investigation, different from that shown in |5]. The main les- 
son that one can infer from these is that boundary conditions 
displayed in Sec. 2 can be significantly weakened. It is clear 
from considerations around <12l l that the point where one can 
expect difficulties is when the adiabatic index F = 5/3, since 
in this case the relevant bounds are relatively poor But the 
numerical data, that are displayed in the forthcoming Table 
and figures, demonstrate suprisingly good agreement with the 
statements of the preceding section. 

TABLE I: Characteristics of the sonic point al, \U,\, R, for the F « 
5/3 poly trope and = 0. 1. The last column shows the areal radius 
of the apparent horizon. The mass accretion rate (second column) 
achieves maximum at mgH ~ 0.663. 



nif/m 




ih 




R. U. al 


Rah 


4.18 X 10- 


-32 


-2.75 X 10- 


48 


2.76 0.424584 0.397394 


1.99 


4.18 X 10" 


-22 


-2.75 X 10" 


38 


2.76 0.424584 0.397394 


1.99 


4.18 X 10" 


-12 


-2.75 X 10" 


28 


2.76 0.424584 0.397394 


1.99 


0.038 




-2.53 X 10" 


18 


2.64 0.424584 0.397395 


1.91 


0.075 




-4.62 X 10" 


18 


2.53 0.424621 0.397703 


1.82 


0.12 




-6.31 X 10" 


18 


2.41 0.424621 0.397702 


1.74 


0.16 




-7.63 X 10" 


18 


2.30 0.424575 0.397359 


1.66 


0.2 




-8.60 X 10" 


18 


2.18 0.424621 0.397702 


1.57 


0.25 




-9.26 X 10" 


18 


2.06 0.424621 0.397703 


1.49 


0.29 




-9.63 X 10" 


18 


1.95 0.424621 0.397719 


1.40 


0.33 




-9.73 X 10" 


18 


1.83 0.424619 0.397784 


1.32 


0.38 




-9.62 X 10" 


18 


1.72 0.424585 0.397395 


1.24 


0.41 




-9.30 X 10" 


18 


1.60 0.424611 0.397586 


1.16 


0.83 




-1.45 X 10" 


18 


0.44 0.424621 0.397707 


0.32 


0.88 




-8.37 X 10" 


19 


0.33 0.424621 0.397707 


0.24 


0.94 




-2.05 X 10" 


19 


0.16 0.424621 0.397704 0.11 



It appears that the characteristics al, Ul, nit/Rt of the sonic 
point practically do not depend, for a given F and a^, on the 
energy density (exemplary results are shown in Table 1), 



suggesting a precision of 10"^. The more reliable measure of 
the numerical error is the difference between the mass to = 1 
(assumed in the equations) and the mass found numerically. 
Our results (see Table ) suggest that numerical error is smaller 
than 0.5 %. The quantity c, is negligible in comparison to 
other sonic point parameters. The mass to, behaves like to, = 
1 - jgaa (Fig. 1). And finally, the extremum of to is achieved 
when to, =s 2/3 (Fig. 2). Notably, numerical results extend 
the regime in which the universality is observed into fluids 
characterized by the speed of sound aco exceeding 1 . This is 
remarkable, since that implies < Rah and the theoretical 
analysis presented above would not be valid. 
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FIG. 1: The dependence of ^ ~ 1 - m, on the asymptotic mass 
density. 
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1 — m» 

FIG. 2: The dependence of m on the ratio (m - m,) /m. 

Formula il6i implies that the mass accretion rate in the case 
with backreaction (when to, < to) can be only smaller from 
the mass accretion rate without backreaction (when to, = to), 
with the same asymptotic parameters poo, a^o, to. Numerical 
results agree with the analytic conclusion that the proportion- 
ality coefficient y depends essentially only on the size of sup- 
port. 
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STABILITY OF SELFGRAVITATING FLUIDS IN 
LAGRANGIAN COORDINATES. 

The effect of selfgravitation of the steadily accreting gases 
on the stabiUty has not been studied, up to our knowledge, 
even in the Newtonian case. In this section we shall fill this 
gap by analyzing nonrelativistic accretion. The relativistic 
accretion will be done elsewhere. The forthcoming analysis 
holds true for polytropes and, more generally, for barotropes. 
The newtonian equations for selfgravitating fluids in terms of 
the Eulerian time t and the areal radius R read 



dtU+UduU = 
dtg = 



Q 



m(R) 



-Ld,{R'eu). 



(18) 



The stability analysis of the Bondi accretion in the Lagrange 
formulation has been performed by Balazs 1,1 31 His approach 
is correct but inconclusive, because of a too stringent under- 
standing of the notion of linearized stability. We show that 
using lagrangian variables one exactly reproduces the stabil- 
ity result already known from the Eulerian approach 01 ■ 

As the primary quantity one takes the deviation ^(f, r) = 
AR{r, t) from the particle position in the unperturbed flow. 
The perturbed velocity reads Af/ = d'f^ = (dt + UdR)(; here 
d'f = dt + UOr is the derivative with respect the Lagrangian 
time. The lagrangian perturbation IS.q can be related to f by 
demanding that Am{R{r - const)) = 0; that implies 



Ag. 



(19) 



The U,a,g solves Eqs. ( I18> with the mass function m(R) 
including the contribution of the fluid itself. One derives from 
the first of Eqs. (I18> the equation ruling ( (take notice during 
the derivation that Am(R(r)) = 0) 



—dup 

Q 



(20) 



The perturbation of the other equation in (I18> is satisfied iden- 
tically. 

Balazs tried to find solutions in the standard way, of the 
form ( (R{r)), t) - exp'"' ({R{r)) with positive co', that is with 
nonincreasing in time modulus; this clearly cannot be done, as 
shown below. The proper approach (adopted in the Eulerian 
variables by Moncrief 1 7]) is to define an energy of perturbed 
modes by the expression 

E ^ J^dVg^ {d,cf + \ {dR(f [a^ - + 



(21) 



where V is an annulus between radii /?, and Roa- The inter- 
esting and nontrivial question (since the last integrand in (" I2U 
can be negative in the vicinity of a sonic point) is whether the 
energy E of M\\ is positive. It is enough to consider the re- 
gion outside of the sonic horizon, because perturbations that 
exist inside will be washed away to the black hole. 
In order to answer that question, notice an identity 

/ 2 -1 9\ ^ r / 2 \ 

g\a - RoRa-j = -5r gRla - ^ I 



2g 



\ 2r) 



■ InR^g; 



(22) 



notice the occurence of the negative term at the end of the for- 
mula. Inserting (I22> into the last term in the integrand of J21> 
and integrating by parts, one obtains after lengthy calculation 
that 



47TRfp{a 



m \ 
2r) 



(23) 



where 



E = dV[x^ + Y^)- 



(24) 



is not definite, while X and Y are given by 

^ = Ve-9r^, (25) 
Y - ^g[-^^f^^.<^-3Rh. (26) 

\R^ Vfl2 - [/2 / 

One can easily show that boundary terms of J23l l are nonneg- 
ative. The occurence of a negative term in i24\ means that the 
energy E of modes perturbing the selfgravitating flow may be 
negative. 

The calculation of the derivative of the energy with respect 
to the Eulerian time yields 

r , -d,m{R) 



R^ 



An pp {d,0RC - t/2) - f/ (5,^)2)]^" . (27) 

A careful analysis of boundary terms shows that the contri- 
bution from Roo and is negative definite. We can conclude 
that the outflow through the inner and outer boundaries is non- 
positive and the energy E of perturbations cannot grow for the 
critical flow. Despite this fact, there is a possibility for unsta- 
ble behaviour of the solutions of the time-dependent hydrody- 
namic equations, because the energy E is not positive definite. 

The situation for subsonic flows is unclear even for test 
fluids, as pointed by Garlick 1 15]. While the analysis at Roc 
yields the same result, that the outflow diminishes the energy 
E, one has to assume that the inner outflow is completely ab- 
sorbed by the surface of an attracting body, that is ^{R\, t) - 0, 
in order to conclude that dtE < 0. For more realistic cases. 



6 



with the occurence of a reflection from the inner boundary, 
the energy of the perturbation may grow, signafling an insta- 
bihty, in addition to the "bulk-type instabihty" suggested by 
the above analysis. 

It is necessary to say in this place that the absence of expo- 
nentially growing linear modes does not mean that perturbed 
solutions of the nonlinear equation will be close to the back- 
ground solution at all times. Linearized stability of a solution 
would imply only that evolving perturbations can be bounded 
by initial perturbations in a suitable sense. The existence of 
exponential linear modes in turn does not necessarily yields 
that perturbations grow indefinitely. Linear instability means 
that some evolving nonlinear perturbations are not bounded 
by initial perturbations; the "strength" of evolving perturba- 
tions does not depend on the "strength" of initial pertrbations 

El. 

Restricting our attention to test fluids, notice the following. 
The energy E is equal to E, and E differs from the expres- 
sion \2Al only by the absence of the last, nonpositive integral, 
hence it is clearly positive. The time propagation in turn is 
Uke in (|23, but without the integral ^^dh^ = ^^-^ 
the right hand side. Thus E cannot grow, which excludes the 
exponential growth of X and Y and the long-term exponential 
growth of linear modes ^. That does not guarantee that \( (R) \ 
is time-independent. 
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